application of indicator and ordinary kringing for modeling of groundwater chloride
نویسندگان
چکیده
introduction the main problems relating to water quality for agriculture are salinity, soil infiltration and specific ions toxicity. accumulation of the specific ions from the irrigation water in the plants reduces crop yields. one of the most common specific ions toxicities results from high concentration of chloride (cl) ion. knowledge of the spatial distribution of cl concentration in groundwater is needed for a better management of the groundwater resources. as limited number of sample data is often available, some appropriate interpolation methods are needed to interpolate between the sample points. ordinary kringing (ok) is a geostatistical estimation method which uses a semivariogram model to predict the unknown values. however, it cannot predict properly the spatial distribution pattern of a highly skewed data. besides, ok estimation variance is not a perfect measure of local uncertainty because it only depends on data configuration not data values. unlike ok, indicator kriging (ik) is a distribution-free approach, which has the ability to model local uncertainty of the estimated values through estimating a conditional cumulative distribution function (ccdf) corresponding to each point. although this method has been used by many researchers for mapping and modeling local uncertainty of various environmental variables such as groundwater cl and other groundwater quality parameters, in iran it has not been employed for such purposes much. the objective of this study is, therefore, to model the local uncertainty of groundwater chloride over kerman plain using ik. the performance of ik is compared with the traditionally used ok (with and without data logarithms). materials & methods study area and sample data this study is performed in kerman plain. the study area is located in a semiarid and arid region. its average elevation above sea level is 1755 m. because of the lack of surface water resources, groundwater resources are the main water resources for agricultural purposes in this area. due to the importance of specific ions toxicity, the groundwater samples were collected from 76 agricultural wells and cl concentration were measured in laboratory. geostatistical analysis first of all, experimental (indicator) semivariograms are calculated to investigate the spatial variability of (within class) cl data. a suitable theoretical model is then fitted to the experimental values for kringing modeling of cl. then ik is used to map groundwater cl and to evaluate the uncertainty attached to the estimates. the results were compared with those achieved from ok and log-kringing (lok). the probability maps of not exceeding two threshold values 10 and 20 meq/lit were generated for cl by ik. these two threshold values are selected according to the irrigation water quality standard proposed by ayres and westcot in 1989. in the following the geostatistical tools and methods used in this study are briefly described: semivariogram the semivariogram quantifies the dissimilarity between observed values as the separation distance between the sample points increases. in practice, experimental semivariogram, γ*(h), is computed for two values separated by a lag distance h as following: (1) where, n is the total number of data pairs of observations z(xi) and z(xi+h)separated by a distance h for a specific direction. kringing needs the semivariogram values for any given lag, therefore, a theoretical model may be fitted to the experimental values and the characteristics of this model can be used. ordinary kringing in ok, the values at unsampled locations are determined by a linear weighted moving averaging of values at sampled locations as: with (2) where z*(x0) is the estimated value of variable of interest at unsampled location x0, li is the weight assigned to the known value of variable at location xi determined based on a semivariogram model and n is the number of neighboring observations. ok produces an estimation variance attached to every estimate, which can be used to generate a confidence interval for each estimate assuming a normal distribution of errors. ok performed on lognormal transformed data is called lognormal ordinary kringing (lok). the estimates have to be back-transformed to the original space at the end. indicator kringing indicator kringing (ik) is based on the coding of the random function z(x) into a set of k indicator random function i(x, zk) corresponding to different cutoffs zk: (3) after transforming the observed data to a new set of indicator variables, the experimental semivariogram is calculated for every set of indicators at each cutoff zk. the conditional cumulative distribution function (ccdf) at each unsampled location, e.g. x0, is then obtained by the ik estimator: (4) where i*(x0; zk) is the estimated indicator transform at unsampled location x0 and li are the weights assigned to the indicator transform i at location xi. these discrete probability functions must be interpolated within each class and extrapolated beyond the minimum and maximum values to provide a continuous ccdf, which covers all possible range of the variable. e-type estimates, which are comparable with ok estimates, may be computed through post processing ik-based ccdfs. local uncertainty measures, e.g. conditional variance and probability maps are also produced through post processing of ik-based ccdfs. in this study kringing methods are performed using the software package gslib. evaluation of the results a cross-validation technique with comparison criteria mean error (me), mean absolute error (mae) and root mean square error (rmse) is used to evaluate the performance of the methods. the most accurate method is the one with the smallest amount of mae and rmse and with a mbe close to zero. results and discussion statistical analysis shows that cl data distribution is strongly positively skewed. a logarithm transform is used to provide a normal frequency distribution of data. experimental semivariograms are calculated for the raw and log-transformed data. for ik, nine thresholds 1.3, 2.4, 3.4, 4.4, 5, 7.4, 12.7, 18 and 24 (meq/lit) corresponding to 10, 20, 30, 40, 50, 60, 70, 80 and 90 percent of cl cumulative frequency distribution functions are selected. then, the observation of cl data is coded according to these selected thresholds and the indicator semivariograms were computed afterwards. the results of semivariogram analysis show that chloride data values and its logarithms are strongly correlated in space and the best fitted semivariogram model is spherical with a range of influence of 42 and 72 km, respectively. however, according to the results, a higher spatial correlation is seen for log-transformed data. cl data had a moderate to strong within class spatial correlation and indicator semivariograms often follow a spherical structure. furthermore, the cross-validation results indicate that lok with the smaller amounts of rmse and mae and a higher amount of correlation coefficient, r, are more accurate than ik for estimating groundwater chloride. beside estimation, one of the main aims of this study is to model local uncertainty of cl data over the study area. both ok and lok provide the uncertainty attached to each cl estimate by calculating its estimation variance. thus, where the estimation variance or standard deviation is smaller, the estimated value of cl is more certain. the results show that ok and lok estimation variance is more related to sampling configuration not to the actual values. in contrast ik conditional variance shows some relation with the sample data in addition to sampling location. besides, ik conditional variance was more appropriate for representing the estimation error than the ok (and lok) variance. moreover, the produced probability maps showed that the probability of chloride exceeding critical thresholds 10 and 20 meq/lit is higher in the northwest and west of the study area. conclusion in this study non-linear indicator kringing is used to model the local uncertainty attached to cl concentration estimates. ordinariy (log) kringing is used to map the spatial distribution of cl estimates. the results show that ordinary log kringing, which is faster and mathematically simpler than indicator kringing, provide more accurate results of cl estimates. the correlation between indicator kringing conditional variance and estimation error is stronger than the correlation between ordinary (log) kringing variance and estimation error. this means that ik conditional ik produces the probability map of not exceeding a critical threshold for cl concentration. these maps can be useful in many decision-making processes, e.g. water resource management.
منابع مشابه
effect of postharvest application of chitosan and calcium chloride on decay and quality attributes of sweet cherry
چکیده ندارد.
15 صفحه اولapplication and construction of carbon paste modified electrodes developed for determination of metal ions in some real samples
ساخت الکترودهاِی اصلاح شده ِیکِی از چالشهاِی همِیشگِی در دانش شیمِی بوِیژه شیمِی تجزیه مِی باشد ،که با در نظر گرفتن سادگِی ساخت، کاربردی بودن و ارزان بودن روش مِی توان به باارزش بودن چنِین سنسورهاِی پِی برد.آنچه که در ادامه آورده شده به ساخت و کاربرد الکترودهاِی اصلاح شده با استفاده از نانو ذرات در اندازه گِیرِی ولتامترِی آهن وکادمِیم اشاره دارد. کار اول اختصاص دارد به ساخت الکترود خمِیر کربن اصلاح شده با لِیگاند داِ...
15 صفحه اولapplication of upfc based on svpwm for power quality improvement
در سالهای اخیر،اختلالات کیفیت توان مهمترین موضوع می باشد که محققان زیادی را برای پیدا کردن راه حلی برای حل آن علاقه مند ساخته است.امروزه کیفیت توان در سیستم قدرت برای مراکز صنعتی،تجاری وکاربردهای بیمارستانی مسئله مهمی می باشد.مشکل ولتاژمثل شرایط افت ولتاژواضافه جریان ناشی از اتصال کوتاه مدار یا وقوع خطا در سیستم بیشتر مورد توجه می باشد. برای مطالعه افت ولتاژ واضافه جریان،محققان زیادی کار کرده ...
15 صفحه اولan application of fuzzy logic for car insurance underwriting
در ایران بیمه خودرو سهم بزرگی در صنعت بیمه دارد. تعیین حق بیمه مناسب و عادلانه نیازمند طبقه بندی خریداران بیمه نامه براساس خطرات احتمالی آنها است. عوامل ریسکی فراوانی می تواند بر این قیمت گذاری تاثیر بگذارد. طبقه بندی و تعیین میزان تاثیر گذاری هر عامل ریسکی بر قیمت گذاری بیمه خودرو پیچیدگی خاصی دارد. در این پایان نامه سعی در ارائه راهی جدید برای طبقه بندی عوامل ریسکی با استفاده از اصول و روش ها...
part a: application of n-(p-toluenesulfonyl) imidazole (tsim) and triphenylphosphine/carbon tetrachloride in several organic transformations part b: application of 8-bromocaffeine for synthesis of some novel 8-caffeinyl derivatives
بخش اول این پایان نامه به طور عمده بر توسعه کاربردهای جدید n-(پارا-تولوئن سولفونیل) ایمیدازول (tsim) و تری فنیل فسفین/ تتراکلرید کربن در تبدیل گروههای عاملی به یکدیگر استوار است. نظر به تنوع زیاد، دردسترس بودن، سمیت کمتر و نقل و انتقال آسان تر الکل ها نسبت به آلکیل هالیدها، تبدیل مستقیم گروه هیدروکسیل به گروه های عاملی دیگر مثل آزید، نیتریل و استر یکی از مهم ترین تبدیلات در سنتزهای آلی است. با ...
15 صفحه اولthe innovation of a statistical model to estimate dependable rainfall (dr) and develop it for determination and classification of drought and wet years of iran
آب حاصل از بارش منبع تأمین نیازهای بی شمار جانداران به ویژه انسان است و هرگونه کاهش در کم و کیف آن مستقیماً حیات موجودات زنده را تحت تأثیر منفی قرار می دهد. نوسان سال به سال بارش از ویژگی های اساسی و بسیار مهم بارش های سالانه ایران محسوب می شود که آثار زیان بار آن در تمام عرصه های اقتصادی، اجتماعی و حتی سیاسی- امنیتی به نحوی منعکس می شود. چون میزان آب ناشی از بارش یکی از مولفه های اصلی برنامه ...
15 صفحه اولمنابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
محیط شناسیجلد ۴۰، شماره ۳، صفحات ۷۵۱-۷۶۴
کلمات کلیدی
میزبانی شده توسط پلتفرم ابری doprax.com
copyright © 2015-2023